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1 Introduction 



Two methods are mainly used at LEP2 to extract the W mass: the thresh- 
old method and the direct reconstruction technique Q. In the first case the 
total W + W~ cross section is measured near threshold (161 GeV), where the 
sensitivity to Mw is stronger, and plotted as a function of the W mass. Con- 
versely, the direct reconstruction method is applied at higher energies and 
requires two steps: 

1. From the experimental data the invariant mass distribution -fe is recon- 
structed. To improve the mass resolution, a constrained fit is usually per- 
formed event by event, assuming no initial-state radiation (ISR) and equality 
between the invariant masses coming from different Ws. 

2. The experimental distribution is compared with the theoretical prediction 
for -jjj and, after Monte Carlo corrections, a reconstructed W mass Mr is 
extracted, with an error AMr. 

Recently, a new method has been proposed (direct fit method), in 
which only the best measured quantities are used to extract the W mass. The 
idea is simple. Given a set of well measured quantities {$} one computes, 
event by event, the theoretical probability V% of getting the observed set of 
values {$j} for {$}. Since this is a function of Mw, given N observed events, 
the logarithm of the likelihood function (L = Hi Pi) is distributed, for large 
N, as a quadratic function of Mw A parabola can then be fitted, from which 
the reconstructed W mass Mr is obtained with an error AMr. 

Although one can always consider more sets {$}, the following choices 
seem reasonable, in practice, for different four-fermion final states |2|, |3|: 

Semileptonic case: i^v^q^q^ 

1. {$} = {E 3 ,fl 3 ,fl 5 ,fl e , Eh}, where Eh is the total energy of the 
two jets. 

2. {$} = {^3,^3,^5,^6}. 

Purely hadronic case: q^q^q^ 

3. {^} = {n 3 ,n 4 ,n,,n 6 }. 
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Purely leptonic case: t^v^ls 

4. {^} = {E 3 ,Q 3 ,E 6 ,Q 6 }. 

In this paper, we explicitly give all the formulae needed to compute the 
probabilities referring to the above four cases and present two FORTRAN pro- 
grams (WEXTER and ERAFITTER) to extract M R from the LEP2 data using 
the direct fit method. The two programs have been developed in parallel 
and continuously cross checked. For this reason we chose to present them 
in a common paper. Also, we decided to put emphasis on the description of 
the common algorithms and on a general illustration of the codes, skipping 
the most technical details. Further information is available directly from the 
authors. 



2 Theory 

Let us start from the case when the analysed events belong to only one 
particular class of four-fermion final states (for example \i~ ~v '^1^2) ■ 

The main problem is computing, for each given event i, the probability 
V({&i}, M w ) of measuring the observed values {$j} for the set {<£>}. This is 
a function of Mw and reads 

where y and a tot are the differential and the total (namely integrated 
over the whole fiducial volume) cross sections for the process under study, 
respectively. Once a to t is known, the computation of y for the four cases 
listed in the introduction is needed. Then, a likelihood function may be 
constructed as follows: 

L = X\V{{^},M W ). (2) 

i 

For large numbers of events, L is distributed as a Gaussian. Therefore, a 
quadratic fit of the form Y = aX 2 + bX + c, with Y = logL and X = Mw, 
can be performed. The maximum of the fitted parabola (— ^-) gives Mr, 

while the statistical error given by the set of analysed events is J—j^- 
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The above procedure can be easily generalized to the case when the anal- 
ysed events belong to different processes. In principle, this allows us to ex- 
tract Mr, using together all LEP2 events []. If the events refer to m different 
processes, the likelihood function is given by 



L = l[V^({<l> l },M ] 



w i 



h3 



P » ({ * a ,M„) = -L-fg, (3) 

where j — 1 : m labels the m different processes. 

Therefore, the most general formula for log L reads 

l°g^ = E Elog (ikx - N ^g^il , (4) 



j=i \i=i 



d{®i}, 



where Nj is the number of analysed events for the j process. 

In the following, we list the algorithms needed to compute ^ or ^ ne 
four cases given in the introduction. Due to the ISR, the sum of the two 
incoming momenta reads 

fs fs 
P = Pi +P2 = (E,p,0,0) , with E = — (xi + x 2 ), p= — {x x - x 2 ) , (5) 

where x\ t 2 are the fractions of energy left to after QED radiation, and 

(pi +P2) 2 = s = xxx 2 s. (6) 



Purely hadronic case 

This is the simplest case because {$} = {f^, Q4, Qq} is 8-dimensional. 
Then, when neglecting ISR, the kinematics of the event is completely deter- 
mined. The inclusion of ISR would instead imply two integrations over x\ 
and x 2 . 

The algorithm is as follows. First one generates x\ and x 2 between and 
1, and computes E and p in eq. (|5|). Then, after parametrizing the final-state 

1 Notice that, since the differential cross sections differ, also CP-conjugate processes 
such as pTv^qiCfr and n + v^qiq2 must be considered separately. 
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momenta as 



Pi = Ei (1, c 0l , s di c^, s di s^) 

i — 3 : 6 , c$ i = cos 6^ , = sin 6>j , etc. , 



(7) 



the four unknown energies Ei are found by using the energy-momentum 
conservation constraints: 



A- 



[ E 3 \ 
V E 6 J 



( 



A = 



1 

C0 3 
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V J 

1 

ce 4 



1 

ce 5 



1 

C0 6 



se 3 C0 3 se 4 c </-4 s o 5 c <f>5 s de,Cct>6 
\ s<? 3 s </> 3 s e 4 s 4 s e 5 s 4>5 s e 6 s<p 6 J 



(8) 



For some values of x\ and x 2 , the above system may give unphysical negative 
Ei. Such configurations must off course be discarded. 

The matrix element squared |M| 2 can then be computed by using the 
reconstructed momenta, so that the kernel cross section, to be convoluted 
with the ISR structure functions, is 



da(x 1 ,x 2 ) 



M 

16 



5 4 (P -P3-P4-P5- Pe) II E i dE i ' 



(9) 



i=3 



where (2ir) 8 and the flux factor have been included in the definition of |M| 2 . 
Since 



J dE 3 dE A dE 5 dE & 5 4 (P - p 3 - p 4 - p 5 - Pe) = , 



the final answer reads 
da(x±, x 2 ) 

dm 



da(xi,x 2 ) E 3 E 4 E 5 Eq 
dVL 3 dVLidVL^dVL^ 16detA 



Ml 



(10) 



(11) 
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Purely leptonic case 

The measured quantities are now {$} = {E 3 ,Q 3 , Eq,Qq}. An additional 
integration over fi 4 is required, besides that over x\ and x 2 . Therefore one 
first generates x±, x 2 , fi 4 and then, from the on-shell condition = (P —p 3 — 
Pi ~Pe) 2 , one gets E 4 : 

= 1 s-2E 3 {E-pc 3 )-2E e {E-pc 6 ) + 2E 3 E e (l- c 36 ) 
4 2 E-pc 4 -E 3 (l- c 34 ) - E 6 (1 - c 46 ) 

Q = C e% , Cjj = COS l(pi,Pj) . (12) 

The four- vector p 5 is given by p 5 = P — p 3 — p 4 — p 6 ; then, for any value of 
x±, x 2 and Q4, the final-state momenta are known (as in the previous case, 
unphysical solutions must be explicitly discarded). The kernel multidifferen- 
tial cross section is therefore 

da(xi,x 2 ) _ da(xi,x 2 ) 



d{<5>i} ' dE 3 dQ 3 dE 6 dn t 



I 



dn 4 — | M)2 _ (13) 

4 16 \E-pc 4 - E 3 (l-c u )- E 6 (1 - c 46 )| ' ' V ; 



Semileptonic case without hadronic energy 

Now {$} = {£"3, fl 3 , fl 5 , Qq} and one integration is needed to compute the 
kernel differential cross section. We chose to integrate over E5. As usual, one 
first generates x\ and x 2 . Then, a bound for E§ is found from the condition 
(P-P3-P 5 ) 2 >0: 

1 s -2E 3 (E -pc 3 ) , , 

< E 5 < ^ — , , • 14 

" 5 ~ 2E-pc 5 -E 3 (l-c 35 ) 1 ' 

By generating £5 in the above interval, computing Eq from the on-shell 
condition (P — p 3 — p 5 — p 6 ) 2 = 0: 

= 1 s-2E 3 (E-pc 3 )-2E 5 (E-pc 5 ) + 2E 3 E 5 (l- c 35 ) ( . 
6 2" E -pc 6 -E 3 (l-c 3e )-E 5 (l-c 56 ) 1 } 

and discarding the unphysical solutions, the final-state momenta are recon- 
structed. Finally, the kernel cross section reads 

da(xi,x 2 ) __ da(xi,x 2 ) 



d{$i} dE 3 dVL 3 dVt b dVL & 
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/ 



= / dE* 



E 3 E 5 E 6 

16 \E - pc 6 - E 3 (1 - c 36 ) - E 5 (1 - c 56 ) 



IMI 



(16) 



Semileptonic case with hadronic energy 

For this case, {$} = {E 3 , fl 3 , Q^, Qq, Eh}, where Eh = E$ + Eq. The set 
{$} is 8-dimensional; therefore, by giving x\ and x 2 , the kinematics of the 
event is completely fixed. The algorithm is as follows. First one generates 
x\ and x 2 , then one computes E 5 and E 6 from the system 



Eh — E 5 + Eq 

(P-P3-P5- P&f 



0. 



(17) 



By imposing E 5:6 — ^ ± 5, one finds the following two solutions: 

S ± = 



2a 

D = /? 2 -4a 7 a = -2(l-c 56 ), 

13 = 2 [p (c 5 - c 6 ) + P 3 (c 36 - c 35 )] , 

7 = s - 2P 3 (E - pc 3 ) - E h (2E -p(c 5 + c 6 )) 



+ £? 3J E? h (2 - c 35 



c 36) + -y (1 - c 56 ) 



(18) 



Then p 4 = P — p 3 — p 5 — p 6 . Not always are both solutions physical. The 
conditions to be fulfilled are 

r 1*1 < f 

D > 
E 4 > . 

Finally, the kernel cross section is 

da(xi,x 2 ) _ da(xi,x 2 ) 



(19) 



d{$i} dE 3 dE h dtt 3 dtt 5 dtt 6 

fd6^E 3 (^-6>)6((P-p 3 -p 5 -p 6 y 

F(5 + ) + F(5~), 
|M(5)| 2 L 



F(S) 



|2a5 + /3| 



(20) 
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3 The program WEXTER 



In this section, we present the first of the two programs that use the described 
probabilities to extract Mr from the LEP2 data. The program WEXTER con- 
sists of three parts: the evaluation of the matrix element, the computation 
of the relevant differential cross sections, by integration over the unobserved 
variables, and the fit to the likelihood curve, to extract Mr. 

The knowledge of the total cross sections a^2 in eq. (P is required as an 
input. The needed a[ J ] can be computed once for all, for each value of M w 
used in the fit, using, for example, EXCALIBUR M. 

After an initialization in SUBROUTINE SETPRO, the matrix element evalu- 
ation is performed in SUBROUTINE MATRIX and SUBROUTINE DIAGA. The dif- 
ferential cross sections are evaluated in SUBROUTINE DIFF, while all needed 
integrations and the fit are performed in the MAIN of the program. 

3.1 The MAIN 

In the MAIN, the input file is read (see later for a detailed discussion). The 
input file must include the name of the data file (DATANAME) containing the 
events to be fitted. In DATANAME, all events must be given in terms of a 
complete set of four-momenta readable with the following format: 

open (unit=2 , f i 1 e =D ATAN AME , status= ' old') 
read (2,60) lp, (p(0,kl) ,p(l,kl) ,p(2,kl) ,p(3,kl) ,kl= 3,6) 
60 format (i4/, (4dl9. 10)) 

where lp is a flag that defines the process and the array p (0 : 3 , 1 : 6) contains 
the four-momenta. 

The first index in p (0:3, 1:6) refers to the component while the second 
one labels the particles: 1 is the incoming e + , 2 the incoming e~, while 3, 4, 
5 and 6 are the four outgoing fermions. 

By convention, the beam is along the x-axis (component 1), with the 
incoming e + along the positive values. 

From the i th event given in the above form, the program reconstructs the 
values for the set {$} of best measured variables (see section 2 and ref. 
[0]), namely energies and angles for charged leptons and solid angles for the 
quarks. The set {$} is automatically determined according to the value of 
LFLAG returned by SUBROUTINE SETPRO (see next section). 
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Of course, the above input format is not suitable when analysing real data: 
in that case the values of the measured variables should be directly 

given as an input. This requires a trivial change in the reading format of 
the program. However, we chose to feed the program directly with the four- 
momenta in order to facilitate Monte Carlo studies. 

Then, the probability "P({$j}, My/) of measuring the observed values 
for the set {$} is computed, by a Monte Carlo integration over the 
initial-state QED radiation - implemented as in ref. || - and, when neces- 
sary, over the unobserved quantities. 

Finally, in the last part of the MAIN, the quadratic fit Y = aX 2 + bX + 
c described in the previous section is performed to extract Mr from the 
analysed events. Also the correlation matrix is computed, to estimate the 
error on the fit due to the Monte Carlo integration. 

3.2 The subroutines SETPRO, MATRIX and DIAGA 

These three subroutines and the whole strategy for the computation of the 
matrix element are taken from EXCALIBUR @|, to which we refer for further 
details. 

In the first part of SUBROUTINE SETPRO the strong and the electroweak 
parameters used in the program are set. They are Mz (ZMI), Tz (WZI), 
sm 2 8 w (STH2), a e . m . (ALPHA) and T w (WWI). Also a s is an input. In the 
program, two different a/s are used. The first one (ALS) controls the coupling 
of the additional gluonic diagrams appearing in four-quark final states || 
(setting ALS= switches off such diagrams). The second one (ALSN) is used 
in the computation of the so-called "naive" QCD factor . 

By default, the widths of the gauge bosons (set in FUNCTION CM2) are 
taken to be fixed and different from zero also in the t-channel. This ensures 
QED gauge invariance of the results ||, but induces a small shift in My/. To 
compensate for it, the variable LSHI is introduced. When LSHI is chosen to 
be 1 in the input file, masses and widths are redefined as described in ref. 
[[]]. If LSHI = such a redefinition is not performed. 

Then, the four-fermion processes, chosen in the input file, are read and 
the corresponding Feynman diagrams built up and printed in the output file. 
In SUBROUTINE SETPRO, the chosen processes are also classified according 
to four categories. The variable that controls the classification is LFLAG: 
LFLAG= I (1= 1 :4) corresponds to the I-th choice for {$} described in the 
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introduction. 

Finally, the right permutation of the four particles in the final state is 
assigned to the variables N3, N4, N5 and N6, for later use in SUBROUTINE 
DIFF. This is relevant for leptonic and semileptonic final states only. In fact, 
all four particles in fully hadronic final states are equivalent. 

In SUBROUTINE MATRIX and DIAGA the matrix element squared is com- 
puted, using the helicity techniques described in ref. Q. 

3.3 The subroutine DIFF 

In SUBROUTINE DIFF the differential cross section is computed according to 
the value of LFLAG given by the event. The algorithms for the computation 
have already been described in section 2. SQJAC is the value returned by 
the subroutine. It is the product of the matrix element times the relevant 
Jacobian. SUBROUTINE DIFF is called from the MAIN of the program, where 
the numerical integration over ISR and not measured variables is performed. 

3.4 The input 

The meaning of the input parameters to be specified in order to run WEXTER 
is the following: 

• OUTPUTNAME (CHARACTER* 15) 

The name of the output file. 

• DATANAME (CHARACTER* 15) 

The name of the file containing the events. 

• NLP (INTEGER) 

The number of different processes contained in the file DATANAME. 

• PAR (3, I) (CHARACTER* 2) 

Produced fermion with label 3 for process I ( to be chosen among 'EL', 
'NE', 'MU', 'NM', 'TA', 'NT', 'DQ', 'UQ', 'SQ\ 'CQ', 'BQ', 'TQ). 

• PAR (4, I) (CHARACTER* 2) 

Produced antifermion with label 4 for process I. 
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• PAR (5, I) (CHARACTER* 2) 

Produced fermion with label 5 for process I. 



• PAR (6, I) (CHARACTER* 2) 

Produced antifermion with label 6 for process I. 

The block of the previous four entries, specifying the processes, should be 
repeated NLP times (1= 1:NLP). 

• NEV (INTEGER) 

The number of events in the file DATANAME to be analysed for the fit. 

• NPF (INTEGER) 

The number of different values of My/ used to fit the likelihood curve. 

• N (INTEGER) 

The total number of integration points. 

• KREL (INTEGER) 

It selects the Feynman diagrams. If KREL = all possible Feynman 
diagrams contributing to the chosen processes are taken into account. 
If KREL= 1, only the CC03 diagrams leading to the reaction e + e~ — > 
W+W~. 

• LQED (INTEGER) 

It includes (LQED= 1) or excludes (LQED = 0) ISR. 

• LCOUL (INTEGER) 

It includes (LC0UL= 1) or excludes (LCOUL = 0) the Coulomb factor 



described in ref. UTO 



LQCD (INTEGER) 

It includes (LQCD= 1) or excludes (LQCD = 0) the "naive" QCD factor 



LENER (INTEGER) 

For semileptonic processes, the sum of the hadronic energies is used 
(LENER= 1) or not used (LENER = 0) as an input. In other words, if 
LENER= 1, a semileptonic process is classified with LFLAG= 1. Other- 
wise with LFLAG= 2. 
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• LFOLD (INTEGER) 

For semileptonic (hadronic) processes a 2-folding (24-folding) over all 
possible jet assignments is performed (LF0LD= 1) or not performed 
(LF0LD= 0). 

• LSHI (INTEGER) 

A shift of masses and widths as described in ref. || is performed (LSHI= 
1) or not performed (LSHI= 0). 

• ROOTS (REAL*8) 

The total energy (in GeV) of the colliding e + and e - . 

• X(I) (REAL*8) 

NPF values of M w (in GeV) used to fit the likelihood curve (1= 1 : NPF). 

• SIG(I.J) (REAL*8) 

NPF values (1= 1:NPF) of the cross section integrated over the fiducial 
volume for the process labelled with J. The I-th cross section must be 
computed with the I-th value for M w (X(I)). These and the following 
quantities are an input for WEXTER and should be computed once for 
all using, for example, EXCALIBUR ||. 

• DSIG(I.J) (REAL*8) 

The NPF values (1= 1:NPF) of the error corresponding to the above 
quantities. 

The block of the previous two entries, specifying the cross sections and their 
errors, should be repeated NLP times (J= 1:NLP). 

3.5 The output 

Since presenting a complete test run output would require the specification of 
too long a list of numbers, such as the events contained in the file DATANAME, 
we decided not to include it here. We just describe what a typical output 
file looks like. After printing out information about the process and the 
parameters given in the input file, the program writes down, for each value 
of M w chosen for the fit, the following quantities: 



j=l i=l 






(21) 
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Then the final result of the fit is reported in the following form: 

FITTED MW WITH A 3 PARAM FIT (Y= a*X"2+b*X+c) : 

wm= 0.803187D+02 +/- 0.503275D-01 +/- 0.750389D-02 

Chi"2/d.o.f. = 0.129242D+01 

a = -0. 197406D+03 
b = 0.317107D+05 
c = -0.128700D+07 

where wm is the reconstructed W mass (Mr) and the first and second errors 
are the statistical and the Monte Carlo errors, respectively. The previous 
result has been obtained by analysing a set of 1600 unweighted CC03 semilep- 
tonic events, produced by EXCALIBUR, with an input mass M w = 80.35 GeV 
at y/s = 190 GeV. This set of unweighted events, together with the input file 
used to get the above output, are available, upon request, from the authors. 



4 The program ERAFITTER 

The program provides all necessary elements to fit the W mass and it is 
based on ERATO ||11|| . Below we give a brief description of the program. 



4.1 The Computational Tree 

The program main . f evaluates the differential cross section for the different 
cases described above. Phase-space generation proceeds through the appro- 
priate algorithms ALG001 to ALG004. Then the momentum assignment is 
used to evaluate the matrix element through the routine MASTER, which is 
extracted from ERATO, and it is specific for each selected channel, i.e. leptonic, 
semileptonic and hadronic. The total cross sections, needed for the evalua- 
tion of the likelihood function, are calculated in the programs evucLmass . f , 
llll_c2_mass . f and qqqq_mass.f. Finally fit2.f collects all information 



and does the actual fit to extract the W mass by using standard MINUIT [12 
calling sequences. 
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In main.f, the input data are read and all physical constants needed for 
the computation are defined. The output is the value of Y^iLi log(da/d{$i}) 
for the various values of M w specified by the user. 

In the subroutines setcl-setc4 the transformation from the Cartesian 
representation of the four momenta to the polar one, needed to define the 
well-measured variables, is performed. 

In the subroutines ALG001-ALG004, the generation of the full phase space, 
including, if specified, the ISR, is performed, according to the discussion 
in section 2. Each subroutine corresponds to the four cases defined in the 
introduction. Then the matrix element is computed by the standard call of 
ERATO. All needed routines are incorporated in the file comrou.f. 

As far as the common blocks used in this code are concerned, they are 
identical to the ones used by ERATO, with the addition of VARI01-VARI04, 
which give the values of the well measured variables in each case, and masses, 
where ndim is the total number of Mw values used in the computation and 
dwmas is an array containing these values. Finally, iaprox defines the ap- 
proximation scheme according to which the Mw dependence is taken into 
account. If iaprox= 1 the Mw dependence is computed by only considering 
the Breit-Wigner functions, 

l/((si 2 - M 2 w f + M W T 2 W ) x l/((s 3 4 - M w ) 2 + M W T W ) , 

which is exact in the CC03 class of Feynman diagrams, and constitutes a good 
approximation for CC10, CC11 and, depending on the cuts, for CC20 as well. 
On the other hand, if iaprox= the exact computation is performed. 

4.2 Input description 

The input needed to run the main . f code looks as follows: 

3 lease l=lb 2=la 3=2 4=3 

190 ! energy 

1 1600 ! first and last generated event read 

0.3 10000 !derr,maxmc 

9 79.73 79.93 80.03 80.13 80.23 80.33 80.43 80.53 80.73 

!nim, dwmas 

1 ! isr 

1 ! ipro icoulomb iaprox 
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' /users/papadopo/Fortr/tmp/newf it/gen/run/evud/unwl90100 . data' 
'lililOOf ; 
'masslOOf ' 

128.07 0.2310309 91.1888 2.4974 80.23 2.033 ! input parameters 
24 !# assignments (foldings) default =1 -> no folding 

The FORTRAN variables corresponding to the above input file are: 

• ial 

The flag defining the choice of the set {$}, according to the order given 
in the introduction. 

• eO 

The collision energy. 

• nevl, nev2 

The first and last event read for analysis. 

• derr, maxmc 

The relative error required in the computation of da/dQ and the max- 
imum number of MC iterations. 

• nim,dwmas(20) 

The number of and the actual values of M w used in the calculation. 

• isr, ipro, icoulomb, iaprox 

The flags for the ISR, the actual process (i.e. e + e~ — > gig2<?3<?4), 
Coulomb correction and the approximation scheme as defined above. 

• pupu 

The input file where the events are stored. 

• pupu 

The output file containing the log-likelihood values for each Mw- 

• pupu 

The output file containing the invariant masses of the W decay prod- 
ucts. 

• ALPHA 1, SINW2 , ZMAS, ZGAMA, WMAS, WGAMA 

The physical constants needed for the computation. 
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• nfol 

The parameter controlling the number of foldings to be done, nf ol=2 
for semileptonic and nfol=24 for hadronic. By folding we mean all 
possible different jet assignments. 

For fit2.f the first input file is: 
inf _min_l : 

' . ./total/qqqqlOO' ! input file (tot_cs or lili) 

'tot2' ! output file (tot_cs_coef ) 

80.23 !mass_0 

! dnorm 

(events 

1,5 !ndatl-ndat2 

0,1 (option 0(2) ,ifirst 1(2) 

with the following correspondence to the FORTRAN variables: 

• FILNAM 

The input file including the total cross section for different values of 
the W mass. 

• FILNAM 

The output file containing the coefficients of the fit of the total cross 
section: 

a = par(l) + par (2) M w + par (3) M w . 

• rmax 

The assumed central value for the measured Mw, usually the input 
value in a Monte Carlo simulation. 

• dnorm 

A normalization parameter (default = 0). 

• nev 

Not used. 

• ndatl, ndat2 

The first and the last rows to be read by the code from the input file. 
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• ioption, ifirst 

ioption=0 and if irst=l in order to perform the total cross section 
fit. 



The second input file reads 



inf _min_2 : 
'lililOOf ; 
'tot2' 
80.23 


1600 
145 153 
2 2 



! input file (tot_cs or lili) 

! output file (tot_cs_coef ) 

!mass_0 

! dnorm 

! events 

!ndatl-ndat2 

Ioption 0(2) , ifirst 1(2) 



with 



• FILNAM 

The input file including the log-likelihood values for each event and for 
the different values of the W mass required. 

• FILNAM 

The output file containing the coefficients of the fit and the fitted M w . 

• rmax 

The assumed central value for the measured Mw- 

• dnorm 

A normalization parameter (default = 0). 

• nev 

The number of events read by the main.f code. 

• ndatl, ndat2 

The first and the last rows to be read by the code from the input file. 

• ioption, ifirst 

ioption=2 and if irst=2 in order to perform the fit to the logarithm 
of the likelihood function. 



17 



The first run performs the fit to the total cross section, whereas the second 
one uses the previous results to extract Mr, including the information from 
the differential cross-section. 



4.3 The output 

Using the input files described above, we performed the fit to a selection 
of 1600 CC03 unweighted events produced by ERATO with ISR and M w = 
80.23 GeV. Then a 24-folding was performed over all possible jet assignments. 
Moreover, the total cross section has been evaluated using the same settings 
(CC03, ISR). The result of the MINUIT program is as follows: 

** 3 **MIGRAD 



MIGRAD MINIMIZATION HAS CONVERGED. 



FCN=.3168413E-01 FROM MIGRAD STATUS=CONVERGED 

13 CALLS 137 TOTAL 

EDM= .57E-21 STRATEGY=1 ERROR MATRIX UNCERTAINTY . 0% 



EXT PARAMETER STEP FIRST 

NO. NAME VALUE ERROR SIZE DERIVATIVE 

1 AO 2278.3 2.8190 -.10390E-09 . 19360E-04 

2 Al 23.469 7.0501 . 16109E-09 .25818E-05 

3 A2 -503.17 22.433 .11872E-08 .25067E-05 



In this output AO, Al and A2 are the coefficients of the quadratic form Y = 
A2X 2 + A1X + AO, with Y = logL and X = (M R - 80.23), from which one 
can obtain Mr in GeV, 

M R - 80.23 = 2.33E - 02 ± 3.15E - 02 ± 7.08E - 03 , 

where the first is the statistical and the second the Monte Carlo error respec- 
tively. 
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5 Conclusions 



We have introduced the relevant formalism and described two FORTRAN pro- 
grams to reconstruct the W mass at LEP2 using the direct fit method intro- 
duced in ref. [0. 
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